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Abstract 


The daytime cloud fraction derived by the Clouds and the Earth’s Radiant Energy Sys- 
tem (CERES) cloud algorithm using Moderate Resolution Imaging Spectroradiometer 
(MODIS) radiances over the Arctic from March 2000 through February 2004 increases 
at a rate of 0.047 ± 0.041 per decade. The corresponding top-of-atmosphere (TOA) short- 
wave irradiances derived from CERES radiance measurements show no significant trend 
during this period. These results suggest that the influence of reduced Arctic sea ice cover 
on TOA shortwave radiation is reduced by the presence of clouds and compensated by 
the increase in cloud cover. The cloud fraction and TOA shortwave irradiance over the 
Antarctic show no significant trend during the same period. 

1. Introduction 

Passive microwave measurements from satellites indicate that the Arctic perennial sea ice 
extent decreased at a rate of 6.4% per decade from 1978 through 2000 [Comiso 2002], 
A more recent analysis by Stroeve et al. [2005] reveals that the September ice extent is 
decreasing even faster, at a rate of 7.7% per decade. Rothrock et al. [1999] show that 
the mean draft of Arctic sea ice has decreased by about 40% between 1958-1976 and 
1990s. All these studies indicate that Arctic sea ice cover is decreasing rapidly. Because 
a changing sea ice cover affects surface albedo, understanding the impact of this change 
on the top-of-atmosphere (TOA) shortwave and longwave irradiances and other related 
atmospheric properties is important for assessing changes in the Arctic radiation budget. 
The Clouds and the Earth’s Radiant Energy System (CERES) instruments [Wielicki et 
al. 1996] and Moderate Resolution Imaging Spectroradiometer [MODIS, King et al. 1992] 
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on NASA’s Terra and Aqua platforms have obtained measurements since March 2000 and 
August 2002, respectively. These instruments provide an opportunity to investigate the 
TOA irradiance and atmospheric property changes associated with the recent period of 
rapid change in sea ice cover. 

Identifying clouds over bright sea ice surfaces and estimating irradiances from radiance 
measurements over polar regions are challenging obstacles to monitoring TOA irradiance 
and atmospheric property changes. Cloud identification over sea ice, however, is expected 
to have improved with the increased number of spectral radiances on MODIS and the 
capability to detect cloud trends has been enhanced by reliable calibrations of all MODIS, 
and CERES channels [Loeb et al. 2006a; Minnis et al. 2006]. Hollinger et al. [1990] 
have also shown that the Special Sensor Microwave/Imager (SSM/I) is a stable microwave 
system. 

The error in estimating TOA irradiances from radiances measured by the CERES 
instruments is smaller than the error in Earth Radiation Budget Experiment [ERBE; 
Barkstrom 1984] irradiance estimates. This is largely due to two factors: better scene 
identification and better angular distribution models. The scene within an approximately 
20 km CERES footprint is identified at a smaller scale (~ 1 km) using collocated MODIS- 
derived cloud and aerosol properties. SSM/I-derived snow and sea ice fractions are also 
collocated within the CERES footprints, which helps in identifying seasonal snow and 
sea ice cover. The new angular distribution models take advantage of improved scene 
identification, resulting in more accurate irradiances than estimated from ERBE [Loeb et 
al. 2006b], Improvements in the irradiance estimates from ERBE to CERES are significant 
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over polar regions [Loeb et al. 2006b], as ERBE had known problems in distinguishing 
clouds from snow and sea ice. In this study, we investigate how TOA shortwave irradiance 
and cloud cover vary with the sea-ice cover changes that have occurred in the Arctic since 
2000 . 

2. Methods 

We use data collected by MODIS and the CERES cross-track instruments on NASA’s Terra 
satellite from March 2000 through February 2004 and on Aqua from August 2002 through 
February 2004. Both Terra and Aqua data are used to determine seasonal variations 
in cloud properties and TOA irradiances, but only Terra data are used for the trend 
analysis. The cloud properties, specifically cloud fraction, in this study are derived from 
MODIS radiances using the algorithms described by Trepte et al. [2002] and Minnis et al. 
[2003]. Hereafter, these products are denoted as cloud products derived from the CERES 
cloud algorithm. The TOA shortwave and longwave irradiances are derived from CERES 
radiance measurements using angular distribution models described in Loeb et al. [2005] 
and Kato and Loeb [2005]. All irradiances and cloud properties are taken from the Single 
Scanner Footprint (SSF) Eddition2BR e vi product for Terra and EdditionlB product for 
Aqua, which are available from the NASA Langley Atmospheric Science Data Center. 

The daily mean TOA shortwave irradiance is estimated from the instantaneous short- 
wave irradiance assuming that atmospheric, cloud, and surface properties are constant 
throughout a 24-hour period. For a given instantaneous irradiance, the corresponding 
daily mean irradiance is given by integrating the irradiance over all daylight hours using 
the solar-zenith-angle dependent albedo from angular distribution models for snow and 
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sea ice from Kato and Loeb [2005] and for land and ocean from Loeb et al. [2003]. No 
correction is made to the longwave irradiance to account for diurnal variations. Changes 
in either shortwave or longwave irradiances over the course of each day due to variations 
in the atmospheric and surface properties are implicitly taken into account because mea- 
surements are taken at many different times of day in polar regions by both Terra and 
Aqua. 

All daily mean irradiances derived from CERES data are averaged over 1° x 1° grid 
boxes and then averaged over a month. To obtain a mean irradiance for 60°N to 90°N, 
hereafter the Arctic, and 60°S to 90°S, hereafter the Antarctic, 1° x 1° regional values are 
area weighted and averaged. To determine the mean cloud fraction over a 1° x 1° grid 
box, the number of MODIS cloudy pixels is divided by the total number of MODIS pixels 
in the box. The cloud fraction over the Arctic and Antarctic is then the area-weighted 
average of cloud fractions over all relevant 1° x 1° grid boxes. Sea ice and snow fractions 
are derived from SSM/I data, which are included in the SSF data set. Similar to the cloud 
fraction, SSM/I-derived snow and sea ice fractions over a CERES footprint are converted 
into an equivalent number of 1-km pixels with snow and sea ice by multiplying the snow 
and ice fraction by the total number of 1-km pixels over the CERES footprint. The snow 
and sea ice fractions over a 1° x 1° grid box are then derived by dividing the number of 
snow and sea ice 1-km pixels by the total number of 1-km pixels in the box. Note that the 
SSM/I-derived snow and sea ice fractions are not available within 50 km of the coastline. 

The monthly mean TOA shortwave and longwave irradiances and cloud amounts for 
the four- year period are derived by averaging the four values from the four years. Anomalies 
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are determined by subtracting the four-year monthly mean value from the monthly value 
for each year. 

3.1. Seasonal variations 

The monthly mean reflected and absorbed shortwave irradiances and longwave irradiances 
show similar seasonal variations for both the Arctic and Antarctic (Figure 1). The month 
with the maximum reflected shortwave irradiance coincides with the month with the max- 
imum TOA downward irradiance for both the Arctic (June) and Antarctic (December). 
This agrees with Curry and Ebert [1992], who show that the maximum reflected short- 
wave irradiance at 80°N latitude occurs in June. It also agrees with the results of Wang 
and Key [2005a], who estimated the shortwave irradiances from the Advanced Very High 
Resolution Radiometer (AVHRR). The maximum (minimum) absorbed shortwave irradi- 
ance occurs in June (December) for the Arctic and December (June) for the Antarctic. 
The absorbed shortwave irradiance in July over the Arctic is nearly equal to that in June 
because the July surface albedo is lower due to melting sea ice. The longwave irradiance 
maximum (minimum) over the Arctic occurs in July (February) and over the Antarctic in 
January (July), approximately one month after the corresponding maxima (minima) in the 
absorbed shortwave irradiance. Net irradiance is positive for two months in the summer 
over both polar regions. 

While the seasonal cycles for the Arctic and Antarctic are similar, these values are dif- 
ferent. The annual average CERES-derived TOA downward shortwave irradiances, short- 
wave albedos, TOA upward longwave irradiances, and TOA net irradiances are shown in 
Table 1, along with the values estimated from the ERBE-likc product ES4 that uses ERBE 
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angular distribution models applied to CERES radiances. Compared to the Antarctic, a 
larger net irradiance during Arctic summer is offset by a smaller net irradiance during 
Arctic winter (Figure lc). As a result, the annual average net irradiances for the two re- 
gions are nearly equal, -92.6 Win -2 for the Arctic and -90.1 Wrrr 2 for the Antarctic. This 
result indicates that the annual average energy transports to the Arctic and Antarctic are 
nearly equal and that about 50% of the energy emitted by the polar regions is supplied by 
dynamical processes that transport energy from mid-latitude regions. 

The CERES cloud fraction retrieved from MODIS shows an increase from approxi- 
mately 50% in winter to about 80% in summer over the Arctic (Figure 2a). The cloud 
fraction maximum occurs in September when the sea ice fraction is lowest. The seasonal 
variation of the Arctic cloud fraction is qualitatively in agreement with the cloud fraction 
derived from AVHRR polar pathfinder data [Wang and Key 2005a] and surface observa- 
tions [Hahn et al. 1995]. The summertime cloud fraction derived by International Satellite 
Cloud Climatology Project (ISCCP, Schiffer and Rossow 1983) and from TOVS [Schweiger 
et al. 2002] over the Arctic for the same four-year period agree with the cloud fraction 
derived by the CERES cloud algorithm. The TOVS-derived cloud fraction is, however, ap- 
proximately 0.3 greater than the cloud fraction from the CERES cloud algorithm in winter. 
This difference could be due to TOVS’ large footprint size (100 km) or greater sensitivity 
to thin clouds. All satellite cloud detections are subject to difficulties in discriminating 
between clear and cloudy conditions in the nocturnal Arctic atmosphere, which is often 
characterized by near-surface temperature inversions. The Arctic winter cloud fractions 
derived from both the CERES cloud algorithm and the MODIS Science Team Collection 4 
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algorithm [King et al. 2003] applied to Terra data, as well as an updated CERES algorithm 
that uses slightly different detection threshold values and channels applied to Aqua data, 
are approximately 40 to 60% and the minimum values occur in January. This contrasts 
with the AVHRR-derived cloud fraction of Wang and Key [2005a], which show a minimum 
during April. 

The seasonal variation of cloud fraction over the Antarctic (Figure 2b) is less pro- 
nounced than that over the Arctic. The mean cloud cover over Antarctica is relatively 
constant, ranging between 62% and 75% during all seasons. The weaker Antarctic sea- 
sonal variation results from somewhat opposite changes over land and ocean. The cloud 
cover over ocean varies from about 70% in winter to about 90% in summer, while for land 
it varies from 65% in winter to 45% in summer. The ocean-land cloud fraction contrast in 
the Antarctic is greater than the contrast in the Arctic, which is not surprising given that 
zonal land-ocean distributions are completely reversed over the two polar regions. 

Beesley and Moritz [1999] suggest that the seasonal cycle of low-level Arctic cloud 
amount is largely due to microphysical processes. Low temperatures during winter allow 
ice crystals to grow larger, resulting in larger fall velocities and a reduction in cloud amount. 
From data taken during the Surface Heat Budget of the Arctic Ocean (SHEBA) project, Lin 
et al. [2003] found that relative humidity increases with surface temperature, especially 
when temperatures exceed 275 K. This implies that the frozen sea surface impedes the 
formation of low-level clouds. Further evidence of this effect is the peak cloud fraction from 
the CERES cloud algorithm occurring in both polar regions coincident with minimum sea 
ice coverage (Figure 2). However, sea ice variability provides only a partial explanation 
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for the land-ocean cloud fraction differences. Other factors that are beyond the scope of 
this paper are also in play and should be explored in the future. 

3.2 Anomalies 

Figure 3 shows time series of the anomalies computed for the Arctic. Because of a large 
change in daytime snow and sea ice fractions, they show a significant (-0.064 ± 0.55 at the 
80% confidence level) trend over the 48 months. The trend in the daytime cloud fraction 
is also statistically significant, but the TOA shortwave irradiance trend is not. Anomalies 
of the same variables examined for the Antarctic do not show significant trends. The 
daytime cloud fraction over the Arctic increases at a rate of 0.047 ± 0.041 per decade 
at an 80% confidence level, which corresponds to a 6.5% relative increase per decade. 
Autocorrelations are taken into account to assess the uncertainty in the trend by reducing 
the number of samples using the method of Wilks [1995]. This rate is larger than that 
derived by Wang and Key [2005b] , who found that the cloud amount derived from AVHRR 
data increased at rate of 0.016 per decade in summer and 0.033 per decade in spring from 
1982 through 1999. However, the latter rates are within the statistical uncertainty of the 
current results and are from different time period. 

The rate of the daytime snow and sea ice cover change derived from this 48-month 
period is -14.4% per decade (the day plus nighttime snow and sea ice cover change is -6.6% 
per decade). If we assume that the difference in the TOA clear-sky albedo between sea 
ice and open water is 0.4 and that the sea ice cover of 45% of the Arctic decreases at a 
rate of 14% per decade in the daytime, the albedo would decrease by 0.024 per decade, if 
clouds are not present. This corresponds to a change in reflected shortwave irradiance of 



-5 Wm -2 per decade. If we assume that the difference of TOA open-water clear-sky and 
cloudy-sky albedo is also 0.1 and clouds do not change the albedo over sea ice, randomly 
overlapping 45% sea ice and 70% cloud covers reduces the TOA shortwave irradiance 
change to -3.9 Wm -2 per decade if the cloud cover stays the same, and to -3.2 Wm~ 2 
per decade, if cloud cover changes at the observed rate. If the clear-sky and cloudy sky 
albedo difference is larger than 0.1, clouds further reduce the irradiance change caused by 
the sea ice cover change. Therefore, the influence of reduced Arctic sea ice cover on the 
TOA shortwave irradiance is reduced by the presence of clouds. The result of the smaller, 
and less statistically significant, reflected shortwave irradiance trend of —2.0 ± 2.0Wm~ 2 
per decade is consistent with the increase in daytime cloud fraction that coincides with 
the steep decrease in sea ice fraction over the Arctic. 

3.3 Ground-based measurements 

To increase confidence in the CERES cloud fractions, we compare them with cloud fractions 
derived from active ground-based instruments. The Atmospheric Radiation Measurement 
(ARM) program operates a millimeter-wave cloud radar and two laser systems, a micro- 
pulse lidar and a Vaisala ceilometer, at Barrow, Alaska. The monthly mean CERES cloud 
amounts derived over the 1° x 1° box that contains Barrow fall between the cloud radar 
and laser derived cloud fractions [Clothiaux et al. 2000] (Figure 4a). The difference of 
cloud radar and laser-derived cloud occurrence is expected because lasers are the most 
sensitive sensors easily available for detecting cloud particles, where cloud radars fail to 
detect smallest cloud particles. The resulting annual mean cloud fractions over Barrow 
retrieved from the cloud radar and lasers are 0.57 and 0.74, respectively, compared to the 
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CERES annual mean cloud amount of 0.71. 


In addition to the ARM observations, NOAA deploys ground-based radiometers to 
measure shortwave and longwave irradiances at Barrow, Alaska, and the Alfred Wegener 
Institute for Polar Research maintains a BSRN site [Ohmura et al. 1998] at NyAlesund, 
Spitsbergen (Norway). Because of the small amounts of precipitable water over polar 
regions, anomalies in cloud fraction should correlate well with anomalies in downward 
longwave irradiance at the surface. Time series of cloud fraction anomalies for the 1° x 1° 
boxes that contain Barrow and NyAlesund are illustrated in Figure 4. The figure also 
shows the corresponding anomalies in the downward longwave irradiance measured at the 
surface. The correlation coefficient of CERES cloud fraction and downward longwave 
irradiance anomalies is 0.57 for Barrow and 0.40 for NnAlesund. 

4. Summary and Conclusions 

Terra data collected from March 2000 through February 2004 suggest that the Arctic 
daytime cloud cover increased by 6.5% per decade during this period while the daytime 
snow and sea ice fraction decreased at a rate of 14.4% per decade. The TOA reflected 
irradiance shows no significant trend during this period, which is consistent with a large 
cloud cover in summer and with increasing cloud cover with diminishing snow and sea 
ice cover. While four years of data are usually not enough to detect trends, the CERES 
data set shows a trend in cloud cover at the 80% confidence level because a relatively large 
change occurs in this period. The large differences in the CERES and TOVS cloud fractions 
suggest that optically thin clouds might often be present during the Arctic winter. This 
result also suggests that the ranges of cloud optical thicknesses retrieved from different 
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instruments needs to be characterized in order to compare cloud fractions from different 
instruments. 

Although based on a limited time series, these results suggest that changes in Arctic 
sea ice are compensated by changes in cloud cover as a result of enhanced evaporation from 
the sea surface, therefore, leaving the TOA energy budget unchanged. The implications are 
that any ice-albedo feedback could be dampened because of increased cloud cover and such 
responses should be sought in climate simulations. However, much additional work needs 
to be conducted to confirm these limited results. To reduce the uncertainty in estimating 
the trend in TOA irradiance and cloud fraction, a longer period of data is necessary. The 
CERES project will be producing an Edition 3 data set beginning in late 2006. It will be 
processed with improved angular distribution models and cloud identification algorithms. 
Thus, further analysis of Arctic cloud cover and radiation trends is left for future research 
when a better and longer data set, especially for nighttime, becomes available. 
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Table 1: Annual Mean TOA Irradiances and Albedos 


TOA Irradiances and Albedo 

Arctic (60°N-90°N) 

Antarctic (60°S and 90°S) 

Downward Shortwave (Wm -2 ) 

204.2 (201.3) 

206.8 (202.5) 

Shortwave Albedo 

0.469 (0.487) 

0.566 (0.591) 

Upward Longwave (Wm -2 ) 

201.1 (201.2) 

179.8 (179.8) 

NET (Wm -2 ) 

-92.6 (-97.9) 

-90.1 (-97.0) 


Numbers in parenthesis are from ERBE angular distribution models applied to CERES 
data. The differences in the CERES and ERBE TOA downward shortwave irradiances are 
due to using spatially and temporally coarse ERBE values in the average. 
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Figure 1 . Monthly mean TOA reflected shortwave irradiances (open circles), absorbed 
shortwave irradiance (solid circles), longwave irradiances (open squares) and downward 
shortwave irradiances (open triangles) for a) the Arctic (between 60°N and 90°N), for 
b) the Antarctic (between 60°S and 90°S), and c) Net irradiance (absorbed shortwave 
irradiance - longwave irradiance) at the top of the atmosphere; the months in parenthesis 
are for the Antarctic. The vertical bars represent the maximum and minimum values for 
a given month as derived from Terra and Aqua CERES data. 
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Figure 2. Monthly mean cloud fractions derived from MODIS radiances by the CERES 
cloud algorithm for a) the Arctic, b) the Antarctic, and c) over ocean surfaces and over 
land surfaces. Monthly mean snow and sea ice fractions derived from SSM/I are shown 
by small solid circles. The cloud amounts estimated from TOVS (Schweiger et al. 2002) 
over the Arctic and International Satellite Cloud Climatology Project (ISCCP, Schiffer and 
Rossow 1983) are shown by closed triangles and open diamond, respectively. The months 
in parenthesis are for the Antarctic and the error bars are the same as in Figure 1. 
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Figure 3 Reflected TOA shortwave irradiance (SW), daytime cloud fraction (Cld. Day), 
and daytime snow and sea ice fraction (Snow Day), anomalies as a function of month for 
the Arctic. Terra data are represented by the solid line with dots, and Aqua data by the 
light solid line with open circles. Open squares indicate months with missing days. The 
horizontal lines are linear regression fit to Terra data. 
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Figure 4 a) Monthly mean day and nighttime cloud fractions over Barrow, AK between 
March 2000 and Feb. 2004. The error bars indicate the maximum and minimum values 
from ground-based active sensors (Clothiaux et al. 2000) during the 4 years. Day and 
nighttime cloud fraction anomalies derived from b) CERES (closed circles) and c) TOVS 
(closed diamond) over a 1° x 1° box including Barrow, AK (71.32°N, 203.39°W) and d) 
from CERES and e) from TOVS over NyAlesund, Spitsbergen (Norway; 79.83°N, 11.93°E). 
The longwave irradiance anomalies derived from surface measurements by pyrgeometers 
and scaled by 0.005 are also shown by open squares. 
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